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By breaking the inversion symmetry of the three-dimensional Dirac metal NasBi, we realize 
topological insulator (TI) phases in the known compound NaBaBi using ab initio calculations. Two 
distinct TI phases emerge: one phase is due to band inversion between the Bi p and Na s bands, 
and the other phase (under pressure) is induced by inversion of the Bi p and Ba d bands. Both 
phases exhibit Dirac-cone-type surface states, but they have opposite spin textures. In the upper 
cone, a left-handed spin texture exists for the s-p inverted phase (similar to a common TI, e.g., 
BLSes), whereas a right-handed spin texture appears for the p-d inverted phase. NaBaBi presents 
a prototype model of a TI that exhibits different spin textures in the same material. In addition, 
NaBaBi may exhibit the introduction of correlation effects to the topological state owing to the 
existence of d states in the band inversion. 

PACS numbers: 73.20.At, 71.20.-b, 71.70.Ej 


I. INTRODUCTION 

In recent years, topological insulators (TIs) have at¬ 
tracted considerable attention owing to their novel fun¬ 
damental physics and potential applications. 1 ^ Because 
of the nontrivial topology of the bulk band structure, 
TIs exhibit gapless surface states inside the bulk energy 
gap. The spin direction is locked to the momentum by 
spin-orbit coupling (SOC). Thus, the topological surface 
state (TSS) exhibits a helical spin texture in momentum 
space. To date, left-handed spin texture is commonly 
found in known TI materials such as the Bi 2 Se 3 fam¬ 
ily 3 and other materials,® 10 whereas right-handed spin 
texture has been reported only in HgS. 11 

Topol ogica l Dirac semimetals were discovered recently 
in NasBi 0 12 and Cd 3 As 2 QMlZIThey also exhibit nontriv¬ 
ial topology in the band structure and can be driven into 
other exotic topological phases, such as TIs and Weyl 
metals. 18 For example, the three-dimensional (3D) Dirac 
point can become massive by breaking a certain crystal 
symmetry, inducing a TI phase. Following this paradigm, 
we report a TI material, NaBaBi, synthesized in a recent 
experimentwhich is equivalent to breaking the inver¬ 
sion symmetry of NasBi by replacing two Na sites with 
one Ba. Furthermore, NaBaBi can be transformed be¬ 
tween two different TI phases via external pressure, one 
with s-p band inversion and left-handed spin texture and 
the other with p-d band inversion and right-handed spin 
texture. 


II. CRYSTALLOGRAPHIC STRUCTURE AND 
CALCULATION METHODS 

NaBaBi was first synthesized via direct reaction of 
the elements Na, Ba, and Bi at high temperature in 
2004.^ As shown in Fig. 1(a), NaBaBi exhibits a hexag¬ 
onal crystal lattice in the P62m space group (No. 189). 



FIG. 1. (color online) (a) Hexagonal crystal lattice 

for NaBaBi. Two atomic layers, [(Bi 3_ ) 2 (Na + )s] 3_ and 
[(Bi 3_ )(Ba 2+ ) 3 ] 3+ , are stacked alternately along the c direc¬ 
tion. (b) Corresponding hexagonal Brillouin zone (BZ) and 
two-dimensional (2D) BZ projected onto the (001) surface, 
(c), (d) Local environments around Ba and Na atoms are oc- 
tahedra and tetrahedra, respectively, distorted by Bi atoms. 


It is a polar intermetallic compound: two atomic lay¬ 
ers are stacked along the c direction in each unit cell, 
one layer with [(Bi 3_ ) 2 (Na + )s] 3_ and the other with 
[(Bi 3- )(Ba 2+ ) 3 ] 3+ . We can understand the lattice of 
NaBaBi as a distorted version of that of NasBi. Two 
equivalent Bi sites that are connected by the inversion 
symmetry in NasBi become non-equivalent in NaBaBi. 
The surrounding Ba and Na atoms are in octahedral and 
tetrahedral environments that are distorted by Bi atoms, 
as shown in Figs. 1(c) and (d), respectively. In prin¬ 
ciple, NaBaBi can also be regarded as a low-symmetry 
version of layered honeycomb ternary compounds, which 
were recently found to be strong 20 ! and weak TIsPI Be- 
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FIG. 2. (color online) Electronic properties of NaBaBi. (a) 
Total and partial density of states. (b,c) Energy dispersions 
along high-symmetry lines in lattice momentum space for 
NaBaBi without and with SOC. Sizes of blue, red, and green 
dots in (b) and (c) are proportional to the weight contribu¬ 
tions of Na s, Bi p, and Ba d orbitals, respectively. Fermi 
energy has been shifted to zero. 


cause the number of total valence electrons is eight in a 
formula unit, the closed-shell configuration indicates that 
NaBaBi is possibly an insulator. 

To investigate the band structure of NaBaBi, we per¬ 
formed electronic structure calculations using density 
functional theory (DFT) with the projected augmented 
wave method, 22 -l as implemented in the Vienna ab initio 
Simulation Package! 23 * 24 ^ The exchange-correlation en¬ 
ergy was considered in the generalized gradient approxi¬ 
mation (GGA) level with the Perdew-Burke-Ernzerhof- 
based density functional! 25 The energy cutoff is set to 400 
eV for the plane wave basis. We adopted k-point grids of 
5x5x7 and 7 x 7 x 9 for the lattice optimization and 
electronic structure calculations, respectively. The sur¬ 
face band structures are calculated on a slab model in a 
tight-binding (TB) scheme based on the maximally local¬ 
ized Wannier functions (MLWFs)P 6 which are projected 
from the bulk Bloch wave functions. 


III. RESULTS AND DISCUSSION 

The band structures in Fig. 2 show that bands near 
the Fermi energy are contributed by the Bi p, Na s, 



FIG. 3. (color online) TB energy band structures of the (a) 
(001) surface and (b) corresponding left-handed spin texture 
around the Fermi energy for NaBaBi. Inset shows Dirac 
points resembling local energy dispersion around the Fermi 
energy at the V points. The surface projected momenta are 
indicated in Fig. 1(b). Fermi energy is set to zero. 


and Ba d orbitals. Without SOC, NaBaBi is a direct- 
gap semiconductor with a band gap of 0.23 eV at the T 
point. The highest occupied band is a Bi p z state, and 
the lowest unoccupied band is a Na s state, as shown in 
Fig. 2(b). When the SOC effect is included, the band 
structure around the T point changes dramatically. SOC 
pushes the Bi p z band above the Na s band, producing 
an inverted band structure with an M-shaped valence 
band [Fig. 2(c)]. Although the Z 2 index cannot be di¬ 
rectly calculated from the parity eigenvalues of the Bloch 
wave functions 27 owing to the inversion asymmetry, the 
SOC-induced band inversion clearly exhibits a topologi¬ 
cal phase transition from a trivial insulator to a TI. 

The TI phase is further confirmed by the direct sur¬ 
face state calculations. Considering that band inversion 
occurs between the Na s and Bi p states near the Fermi 
energy, the MLWFs are derived from the atomic-like Na s 
and Bi p orbitals to obtain the TB parameters. Then we 
construct the Hamiltonians of a slab for the (001) surface 
with a thickness of 200 unit cells and obtain the surface 
band structure from direct diagonalization. As shown in 
Fig. 3, the surface exhibits nontrivial Dirac-cone-type 
TSSs inside the bulk band gap. Because of the absence 
of inversion symmetry, the two Dirac cones on the top 
and bottom surfaces of the slab are not degenerate in 
energy, especially when they extend far from the Bril- 
louin zone (BZ) center. As in most TIs, the surface Dirac 
cone of NaBaBi exhibits left-handed spin texture in the 
upper Dirac cone [Fig. 3(b)]. These results unambigu¬ 
ously demonstrate the topologically nontrivial features of 
NaBaBi. 

As observed in the bulk band structure, the Ba 5 d state 
appears as the second-lowest conduction band above the 
Na s state. Next, we show that the compressive pressure 
can shift the Ba d band down to realize an inversion of 
the Ba d and Bi p bands, i.e., a new TI phase. At the 
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same time, the Na s state is pushed up into the conduc¬ 
tion bands, and the original Na s and Bi p inversion is 
removed. In the crystal field of the BaBU octahedron 
[Fig. 1(c)], the low-energy t 2 g bands are further split by 
strong distortion of the octahedron, where Ba d z 2 be¬ 
comes the lowest band of the t^g states (the z direction 
corresponds to the crystal c axis). Under compressive 
pressure that shrinks the BaBU octahedron, the crystal 
field splitting of the Ba d states increases, so the lowest 
band, Ba d z 2 , is pushed down in energy. In the crystal 
field of the NaB4 tetrahedron [Fig. 1(d)], compressive 
pressure will simply increase the energy of the Na s state 
by enhancing the Coulomb repulsion of Bi 3- anions. 

In the calculations, we apply the external pressure in 
two different ways: as hydrostatic pressure and as uni¬ 
axial pressure. Under hydrostatic pressure, we first es¬ 
timate the bulk modulus by fitting the total energy de¬ 
pendence on the volume with Murnaghan’s equation 28 
and obtained Bq = 25 GPa, which awaits validation by 
future experiments. Under uniaxial pressure, a tensile 
strain is applied along the c axis while the a and b axes 
are optimized, which is equivalent to applying in-plane 
compressive strain in the ab plane, such as that induced 
by a substrate. Under both types of pressure, the Ba 
d z 2 band shifts down and the Na s band shifts up, as 
expected. As a consequence, two topological phase tran¬ 
sitions occur with increasing pressure. During the first 
transition, at a hydrostatic pressure of 1.5 GPa (uniaxial 
pressure of 2 GPa), because of the upshift of the Na s 
state, the band inversion between Na s and Bi p disap¬ 
pears. Thus, NaBaBi changes from a TI into a normal 
insulator. During the second transition, at a hydrostatic 
pressure of 7 GPa (uniaxial pressure of 12 GPa), because 
of the downshift of the Ba d xy state, a new inversion 
occurs between Ba d and Bi p at the T point. There¬ 
fore, NaBaBi changes from a normal insulator to a TI 
phase. We show the phase diagram with respect to the 
pressure in Fig. 4, where two TI phases, nontrivial-I and 
nontrivial-II, are separated by a trivial insulator phase. 

As an example, we perform a detailed analysis of the 
electronic structures at a hydrostatic pressure of 10 GPa. 
The nontrivial-II phase at other pressures (hydrostatic 
and uniaxial) exhibit the same band order and just dif¬ 
ferent band gaps. The DFT band structures without 
and with SOC are compared in Figs. 5(a) and (b). Be¬ 
fore SOC is included, the Na s state is pushed up as 
the second-lowest conduction band compared to the case 
without pressure, and Ba d? z becomes the lowest conduc¬ 
tion band. Then SOC induces band inversion between 
the Ba d 2 z and Bi p z states. The evolution of the NaBaBi 
band structures is further illustrated with respect to SOC 
and pressure. To demonstrate the topological feature, we 
also calculate the surface state on the (001) surface us¬ 
ing a slab model. The MLWFs here are derived from the 
Na s, Bi p, and Ba d atomic-like orbitals. The corre¬ 
sponding surface band structure is shown in Fig. 6(a). 
We can clearly see the Dirac-cone-type TSSs inside the 
bulk band gap. Unlike the previous TI phase with s-p 


(a) uniaxial strain 
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FIG. 4. (color online) Evolution of band inversion strengths 
of Es- P = E^a-s — E^i- Pz and E d _ p — E] 3a _ dz2 — E^ i _ Pz 
under (a) uniaxial strain in c direction and (b) hydrostatic 
pressure. Negative values of the inversion strength represent 
inverted band gaps, i.e., the TI phase. Vo and Co are the 
volume and lattice constant in the equilibrium state, respec¬ 
tively, whereas V and c are those under the pressure or strain 
state. The two shaded areas in both top and bottom panels 
are topologically nontrivial states with p-d (nontrivial-II) and 
s-p band inversion (nontrivial-I), and the bright blocks are 
topologically trivial states. 


inversion, however, this p-d-inverted TI phase exhibits 
right-handed spin texture for the upper Dirac cone on 
the top surface, as shown in Fig. 6(b). Therefore, we ob¬ 
tain different spin textures in a single material, NaBaBi, 
via band structure engineering. 

To understand the origin of the different spin tex¬ 
tures in NaBaBi, we also construct the effective Hamil¬ 
tonian of this system using the invariant theoryP 9 We 
focus on the low-energy physics around the T point, at 
which the wavevector group is D^h- Because the spin 



































4 



FIG. 5. (color online) Bulk band structures for NaBaBi (a) 
without and (b) with SOC under a hydrostatic pressure of 10 
GPa. Sizes of blue, red, and green dots represent the weight 
contributions of Na s, Bi p, and Ba d orbitals, respectively. 
High-symmetry momenta are implied in Fig. 1(b), and the 
Fermi energy has been tuned to zero. Original s-p inversion 
changes to p-d inversion. 


should be responsible for the opposite signs of ^ and 
In the appendices, we also discuss the relationship 
between the spin texture and mirror Chern number for 
the s-p and p-d inversions. This relationship is similar to 
that for the opposite spin textures in HgTe and HgS . 11 



is taken into account, we need to consider the double 
group of Dsh, which possesses only three 2D irreducible 
representations (irreps), denoted as r 7 ; 8 ; g, which are de¬ 
fined in Appendix A. It turns out that both the s or¬ 
bital of Na atoms and the d z 2 orbital of Ba atoms be¬ 
long to the f 7 irrep, whereas the p orbital of Bi atoms 
corresponds to the f 8 irrep. This suggests that both 
the s-p inversion and p-d inversion can be described 
by the same effective Hamiltonian. Under the basis 
|H 7 ,1/2), |f 7 , —1/2), |T 8 , 1/2), |r 8 , —1/2), _where the r 7 
band is for the s or d orbital and the T 8 band is for 
the p orbital, the effective Hamiltonian can be written as 

H a = e(k)/ 4X 4+M(k)r5+P a ^r 4 5+Q a (A; x r25-A: 2/ ri5), 

where e a (k) = e a ^ + e a p + e a?2 ^, and M a ( k) = 
M a ,0 + Here a = s or d when the 

r 7 states are given by the s or d orbitals, respectively. Ti 
and Tij are T matrices, and e, M, P, and Q are material- 
dependent parameters. The details of the construction of 
the model are given in Appendix A. 

For comparison with the ab initio calculations, we 
choose an open boundary condition along the z direc¬ 
tion and calculate the surface states and the correspond¬ 
ing spin texture on a finite slab using the above effec¬ 
tive Hamiltonian. We find that the spin texture is right- 
handed for ^ < 0 and left-handed for ^ > 0 , which is 
shown in detail in the appendices. A comparison of this 
with the ab initio calculations suggests ^ > 0 for s-p in¬ 
version and 7 ^- < 0 for p-d inversion. The parameters P a 

W d 

and Q a are related to microscopic wave functions at the 
T point. P a = -a^(a\d z \p z ), and Q a = -b^(a\d x \p x ), 
where a and b are coefficients for the coupling between 
| Pz,&) and |p±, <j) of the p orbitals due to SOC, as shown 
in Appendix B. The different crystal environments for the 
s orbitals in Na atoms and the d orbitals in Ba atoms 


FIG. 6. (color online) (a) TB surface band structure on (001) 
surface with p-d band inversion under hydrostatic pressure 
of 10 GPa. Inset shows local bands around Fermi energy 
at r point, (b) Corresponding surface spin helical Dirac 
cone. Original left-handed spin texture in p-s inversion phase 
changes to right-handed. 


IV. SUMMARY 

In conclusion, inspired by the 3D topological Dirac 
semimetal Na 8 Bi, we theoretically propose another topo¬ 
logically nontrivial state in NaBaBi. Two distinct TI 
states with different surface spin textures are achieved 
via band engineering. The native TI phase, with s-p 
band inversion, exhibits left-handed spin texture in the 
Dirac-type surface states. Under external pressure, the 
other nontrivial state, which has p-d inversion and ex¬ 
hibits right-handed helical spin texture, can be induced. 
This is further explained by the effective Hamiltonians 
constructed, in which s-p and p-d inverted TI phases ex¬ 
hibit opposite mirror Chern numbers and thus opposite 
spin helicity. Because the d orbitals participate in the 
band inversion, it is possible to bring the correlation ef¬ 
fect into the TI phase in this compound. 
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Appendix A: Construction of effective Hamiltonian 
by symmetry 

In this section, we will show how to construct the ef¬ 
fective Hamiltonian on the basis of the crystal symmetry. 
The point group at the T point in momentum space for 
NaBaBi is D 3 ^. The character table for the D 3 h dou¬ 
ble group is shown in Table [l| To distinguish the T ma¬ 
trices that we will use later, we use f to denote differ¬ 
ent irreducible representations (irreps). Only double ir- 
reps can be used to describe electron systems with spin. 
Thus, we consider only the rV ? 8,9 irreps. The basis func¬ 
tions for the T 7 , 8,9 irreps are denoted as IF 7 ) = |j = 
1/2 ,m,j = ±1/2), |f 8 ) = z\j = 1/2 ,mj = ±1/2), and 
|Tg) = |j = 3/2,mj = ±3/2), respectively. 

TABLE I. Character table for the double group of D 3 h- 


(|r 7 ,1/2), |f 7 , -1/2), |f 8 ,1/2), |f 8 , —1/2)) T . The 

TABLE II. Character table for s, p z , and d z 2 
orbitals with spin. 


Orbitals with spin 

2C 3 

253 

Irreps 

st, si 

1 

V3 

f 7 

Pz t, Pz i 

1 

-V3 

f 8 

d z 2 t, d z 2 4 . 

1 

V3 

f 7 


transformation matrices of the symmetry operations in 
D%h group under this basis are 

• Time reversal operation: T = 

diag[—i(T2K, — icr 2 X] = —73 0 icr 2 A; 

• Mirror operation along the z direction: A4 Z = 
diag[ia 3 , -icr 3 \ = r 3 0 icr 3 ; 

• C 3 operation: C 3 = diag[e ~ 1 ^, e 1 ^, , e z 5]; 


D 3 h 

E 

E 

(Jh&CTh 2(73 

2C3 

2S 3 

25 3 

3 C , 2 & 3 C , 2 

S(J y 3(X y 

• C2,x operation along the x direction: C2, x = 

Tt 

1 

1 

1 

1 

1 

1 

1 

1 

1 

diag[—iai,icri\ = —r 3 0 icr 

r 2 

1 

1 

1 

1 

1 

1 

1 

-1 

-1 


f 3 

1 

1 

-1 

1 

1 

-1 

-1 

1 

-1 

• Mirror operation along the y direction: M y = 

f 4 

1 

1 

-1 

1 

1 

-1 

-1 

-1 

1 

diag\ia 2 V (JC 2 \ = t 0 0 icr 2 , 

Fb 

2 

2 

-2 

-1 

-1 

1 

1 

0 

0 


r 6 

2 

2 

2 

-1 

-1 

-1 

-1 

0 

0 

where and are the Pauli matrices in spin and 

r 7 

2 

-2 

0 

1 

-1 


-Vs 

0 

0 

band space, respectively. Here we select a particular 

r 8 

2 

-2 

0 

1 

-1 

-Vs 

Vs 

0 

0 

gauge such that T = 102 K in spin-up and spin-down 

f 9 

2 

-2 

0 

-2 

2 

0 

0 

0 

0 

space, where K represents complex conjugation. We set 


The orbitals of NaBaBi near the Fermi energy are s 
orbital of Na atoms, p orbital of Bi atoms, and d orbital 
of Ba atoms. From a first-principle calculation, we find 
that the main component of the p (d) orbital of Bi (Ba) 
atoms is p z (d z 2 ) orbital. Next, we will determine the 
irreps to which the 8 , p z , and d z 2 orbitals belong. 

To distinguish irreps 1 ^ 8 , 9 , we need to consider 
only the characters under C 3 and S 3 operations. 
C 3 \s(d z 2,p z ),p = e“ i7r / 3 |s(<42,p z ),t). C 3 \s{d z 2,p z ),l 
) = e l 7 r / 3 \s(d z 2,p z ), 1 ). All the characters under C 3 rota¬ 
tion for the 8 , p z , and d z 2 orbitals are x = e _ 27 r / 3 +e 27r / 3 = 
1 . 

For the S 3 operation, S 3 \s(d z 2 )^) = ie~ Z7r / 3 \s(d z 2 ), t), 
and S 3 \s(d z 2 ),\) = — ie l7r / 3 \s(d z 2 ) : |). The character is 
X = ie -27r / 3 — ie l7r / 3 = y/3. Thus, both the 8 and d z 2 
orbitals belong to irrep T 7 . For the p z orbital, however, 
S 3 |p*,t) = ~ie~™ /3 \Pz,t), and 5 3 |p 2 U) = ie in / 3 \p z ,l). 
The character is x = — ie _47r / 3 + ie m / 3 = — \/3. Thus, 
the p z orbital belongs to irrep T 3 - 

As shown in Table [TI| the s and d z 2 orbitals belong to 
irrep F7, whereas the p z orbital belongs to irrep T 3 . 

After identifying the irreps for each orbital near the 
Fermi energy, we use the representation theory to con¬ 
struct the effective model for NaBaBi near the T point 
in momentum space. Basically, we follow the procedure 
in Ref. [30j 

The effective Hamiltonian we will con¬ 
struct has the following basis: £ = 


By using this convention, we make all the parameters real 
in the Hamiltonian that we will derive. 

We start with 5 Dirac T matrices on the basis men¬ 
tioned above, which are defined as 


Ti = Ti < 

r 4 = t 2 ' 


•oi, 

> 1 , 


r 2 = n 0 cr 2 , 

r 5 = 7-3 01 . 


r 3 = ti 1 


>< 73 , 

(Al) 


They satisfy Clifford algebra: {T a ,T & } = 2 5^. The other 
10 T matrices are constructed by T a ^ = [T a ,Tb]/2z and 
are expressed explicitly as 


T ij = [ti 0 a i: ri 0 aj]/2i 




r i4 = [t\ 0 CTj, r 2 0 l ]/2 i = r 3 ( 8 ) 

T i5 = [ti 0 (Ti, r 3 (8) l]/2* = —r 2 <S> o’, 
T 45 = [ t 2 0 1 , t 3 0 l]/2 i = t\ 0 1 , 


(A2) 

(A3) 

(A4) 

(A5) 


where i, j = 1, 2, 3. 

Under the time-reversal transformation matrix, the T 
matrices satisfy 


TYiT- 1 = Ti, 
TTijT _1 = -Tij, 


i = 1,2,3,4,5; 
= 1,2,3, 
Tr i4 T _1 = —T i4 , * = 1,2,3; 

TIVT- 1 = -r i5 , * = 1,2,3; 

rr^r - 1 = -r 4 


- 45- 


(A6) 

(A7) 

(AS) 

(A9) 

(A10) 


Using the Ad^, C^, and M. y transformation matrices, 
we have 
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M z :M Z r 1>2>5 ,12,34,15,25 A4 2 1 — ri, 2 ,5,12,34,15,25, 

M z r 3A ,23,31,14,24,35,45*^^ 1 — — 1^3,4,23,31,14,24,35,45 ? (All) 

^2,x ^ 2 ,^r 2 , 3 ,5,23,14,25,35^,^ = 1^2,3,5,23,14,25,35? 

C2,^1,4,31,12,24,34,15,45^2,* = —1^1,4,31,12,24,34,15,45? (A12) 

My '-MyT 2 , 4 ,5,31,24,25,45-My 1 = 1^2,4,5,31,24,25,45? 

Al y ri ? 3 ? i2,23,14,34,15,35Al 2/ 1 = — Ti ? 3 ; i 2 ,23,14,34,15,35- (A13) 


Under C 3 rotation symmetry, the T matrices satisfy 
Y' e = e* s / 26 , re _2E / 20 , with E = 1 0 cr z . Furthermore, 

= |[E,r^]. Thus, the commutation relations be¬ 
tween E and T are important here, and are listed as fol¬ 
lows: 


[S.rq 

= 2fr 2 , 

[E? T 2 ] — — 2iTi, 

(A14) 

P,r 3 ] 

= [E,r 4 ] = [s,r 5 ] = 0 , 

(A15) 

[s,r 12 ] 

= 0, 

[E,r 34 ] = 0 , 

(Aie) 

[s,r 3 i] 

= —2 iY 23 

? P? r23] = 2 zT 3 i, 

(A17) 

P,r 14 ] 

= 2iT 24 , 

P? r 24 ] = — 2 iTi 4 , 

(Ais) 

P,r 15 ] 

= 2iT 25 , 

P? ^ 25 ] = — 2iTi5, 

(A19) 

[s,r 36 ] 

= 0 , 

P? 1 ^ 45 ] = 0 . 

(A20) 


Using these commutation relations, we can easily solve 


for Y' e , which are expressed as 

ri ( 0 ) = Ti cos 0 0 T 2 sin 0 , 

T' 2 (0) =risin(9-r 2 cos(9, (A 21 ) 

r ' 23 (0) = r 23 cos o 0 r 3 i sin 0, 

r ' 31 (0) = r 3 i cos 0 - r 23 sin 0, (A22) 

r'i 4 ( 0 ) = Ti 4 cos 0 0 T 24 sin 0 , 

T' 24 (0) = ri 4 sin<9 —r 24 cos<9, (A23) 

r i 5 (#) = r i 5 cos 0 0 r 25 sin 0 , 

r 25 (#) = r i 5 sin 0 - r 25 cos 0, (A24) 

r'(0) = r 3 , ri(0) = r 4 , 

r' 5 ( 0 ) = r 5 , r ' 4 = r 34 , ri 2 = r 12 , 

r ^ 5 = r 35 , r ' 5 = r 45 . (A25) 


T matrix pairs [Ti, T 2 ], [T 23 , T 3 i], [ri 4 ,r 24 ], and 

[r 15 , r 25 ] transform as the x,y components of a vector 
under C 3 rotation. The other seven T matrices trans¬ 
form as scalars or pseudo-scalars under C 3 rotation. 


Because Yj <g> Yg = T 3 0 T 4 0 Ts, here we consider 
only irreps Y 3 , f 4 , and Y$ for the effective model. The 


character table for these three irreps is given in Table III 


Under C 3 rotation, we have two rotationally invari¬ 
ant combinations of T matrices and k that preserve 
time-reversal symmetry. They are k' x r' 2 5 — k' y Y' 15 = 
k x Y 2 5 k y Y \5 and k x Yi§ 0 kyY ?25 • How¬ 

ever, under the M y operation, MyY'^My 1 = —T' 15 , and 
A^yT^Al^ — Tr 2 3 , M y k x — k x , and A/lyky — ky 
The effective Hamiltonian needs to preserve the mirror 
symmetry M y and time-reversal symmetry simultane¬ 
ously. Therefore, the qualified term reads k x Y 25 — k y Y 15 . 

TABLE III. Character table of Y matrices and 
polynomials of the momentum k. 


Representation Basis functions T 


r 3 

r 3 

0 


r 35 

- 

f 4 

k z 

- 


r 4 

0 


r 45 

- 

f 5 

{kx ? ky} 

- 


{kx ky ? k x ky 

) + 


(ri,r 2 ) 

0 


(ri 5 ,r 2 5) 

- 


Similarly, we obtain another term, 2k x k y Y 1 + (k 2 — k 2 )Y 2 . 
Thus, the effective Hamiltonian on the new basis, £ = 
(|f 7 , 1 / 2 ), |f 7 , - 1 / 2 ), |f 8 , 1 / 2 ), |f 8 , - 1 / 2 )) T , up to the 
second order of k, is expressed as H = e(k)/ 4x4 + 
Af(k)T 5 + Pk z Y^ 5 0 Q(k x Y 25 — k y Yi$) 0 Q f (2k x k y Y 1 0 
(k 2 —k 2 )Y 2 ). Explicitly, the effective Hamiltonian can be 
written as 


M 7 ( k) 0 


Heff = 


Pk z 


M 7 ( k) -Qk+ + Q'(2k x k y + i(k 2 x - k 2 )) 


h.c. 


M 8 { k) 


Qk _ + Q’(2k x k v - i(k 2 x - k 2 y )) \ 
Pk z 
0 

M 8 ( k) 


(A26) 


where 


e(k) = e 0 0 0 e 2 k 2 , 

M(k) = M 0 + M 0 + M 2 A£ 
M 7 (k) = e(k) + M(k), 

M 8 (k) = e(k) - M(k), 


and e,, Mi, P, P', and Q are real. Because e(k), 


(A27) 
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M( k), 74 x 4 , and T 5 belong to the fi irrep and they are 
all even under time-reversal symmetry, combinations of 
these terms, e(k)/ 4 X 4 and M(k)T 5 , are invariant under 
all the operations of the D^h point group. Thus, we also 
need to include them in the effective Hamiltonian. This 
Hamiltonian is used to explain the spin textures for both 
the s-p and d-p inversions in NaBaBi. 


Appendix B: Spin helicity and mirror Chern number 

We consider a sla b described by the effective Hamil¬ 
tonian in Eq. ( |A26 ) with the open boundary along the 
k z direction. The spin texture is calculated numerically 


with appropriate parameters for qualitative agreement 
with the result obtained by the first-principle methods. 
The parameters are listed in Table IV In Fig. [7| we show 
the Fermi arc and spin texture on the upper Dirac cone of 
the top surface of NaBaBi on the k x k y plane. The Fermi 
arc is a circle, and the directions of the spins are tangen¬ 
tial to the circle. The spin texture is right-handed (left- 
handed) for NaBaBi with d-p inversion ( s-p inversion), 
which is consistent with the results of the first-principle 
calculations. The key point for the opposite spin tex¬ 
tures is found to be the sign of P/Q. For P/Q > 0 , the 
spin texture is left-handed, whereas for P/Q < 0 , it is 
right-handed. 


TABLE IV. Parameters for NaBaBi. For s-p inversion, a = 0.05 eV - A, Q — 0.5 eV • A, and Q' = —0.2 eV • A 2 ; for d-p 
inversion, e± = —0.05 eV - A, Q — —0.5 eV • A, and Q' = 0.2 eV • A 2 . 

e 0 [eV] ei[eV • A 2 ] e 2 [eV • A 2 ] M 0 [eV] Mi[eV • A 2 ] M 2 [eV • A 2 ] P[eV • A] Q[eV • A] Q'[eV • A 2 ] 

0 ± 0.05 0 -0.5 5 5 1 ± 0.5 ± 0.2 


This can be understand from the mirror symme¬ 
try along the y direction. Next, we hope to unveil 
the physics behind the relationship between P/Q and 
the spin texture. Because the system preserves mir¬ 
ror symmetry with respect to the k y = 0 plane, the 
Hamiltonian commutes with the mirror operator M. y : 
M. y H(k x , 0, k^Aiy 1 = H(k x ,0,k z ). Therefore, the 


Hamiltonian and the mirror operator can be diagonal¬ 
ized simultaneously. The mirror operator along the y 
direction is given by Ai y = To ( 8 ) icrz- In the +i mir¬ 
ror parity subspace, the common eigen-wavefunctions are 

^(+*)i = 75 (M>°>°) t and ^(+*)2 = ^(0,0,1,0 T - 

On these two basis functions, the projected Hamiltonian 
reads 


77 +,= 


m 7 P(k z 

P(k z -i^k x ) + Q'k 2 x 


i%k x ) 

Ms 


Q'k 2 x 


(Bl) 


where Mj(k) = e(k) + M( k), and Mg( k) = e(k) — M( k). 
Q' is much smaller than P and Q and is neglected in the 
analysis below. 

The Hamiltonian in the +i mirror parity subspace has 
a nonzero Chern number, which can be calculated from 

n = Ibz €r d • x mA where d = 1 X Mx,dy,dz) 

when we write T7+, as T7+, = e(k) + d • cr. 3 2 Because 
the Chern number for the effective Hamiltonian in 
the —i mirror parity subspace always takes the opposite 
value because of time-reversal symmetry, = — n+,, 
the total Chern number is zero, and only the mirror 
Chern number, defined as um = n+i ~ n - i = n+J^, can be 
nonzero. For P/Q > 0, um = — 1, whereas for P/Q < 0, 
tim = 1. The mirror Chern number is directly related 


to the spin texture of the surface states according to the 
bulk-edge correspondence. Now let us consider that the 
open boundary along the £ direction and the mirror sym¬ 
metry along the y direction still exist for k y = 0. For the 
mirror parity subspace and n+, = 1 , there is a chiral 
edge mode with positive velocity along the k x direction. 
Because the mirror parity + i subspace corresponds to 
spin up along the y direction, this chiral mode should 
have spin up along the y direction. The spin-down state 
along the y direction is in the —i mirror parity subspace 
and thus possesses negative velocity . 11 Therefore, the en¬ 
tire system has a right-handed spin texture for P/Q < 0 . 
When P/Q > 0, the mirror Chern number is —1, and the 
spin texture is left-handed. 
















(b) 

► 


(a) ^ -► 

1* 

f 

K 




* fc. 


if i 
* ' ->• kx ), 


t 


4 - ^ 

s-p inversion 


d-p inversion 


FIG. 7. (Color online). Spin textures on the Fermi arc 
with energy Ef = 0.1 eV (compared with the Dirac point en¬ 
ergy) on the upper Dirac cone of the top surface of NaBaBi. 
Panels (a) and (b) correspond to s-p and d-p inversion in 
NaBaBi, which have left-handed and right-handed spin tex¬ 
tures, respectively, on the upper Dirac cone of the top surface 
of NaBaBi. 


Next, we will explain what determines the opposite 
signs of P/Q for s-p and d-p inversion. Without losing 
generality, we neglect the off-diagonal quadratic term for 
simplicity. We can write the basis functions for irreps T 7 
and Tg in terms of the 8, p z , and d z 2 orbitals. In the +i 
mirror parity subspace, there are two ways to express the 


basis functions for irrep T 7 , ^(|r 7 , l/2) z +i\T 7 , —1/2)^), 
which are ^(s) +i = i^{\s, t *) + i | s 4 *)) and ^(d z 2 ) +i = 
i^(|c4 2 ,tz) + i\d z *,\- z )), where the spins are along the 
z direction. The basis function for irrep Tg in the mirror 
parity A = +i subspace is 'i>{p z )+i = ^(|r 8 ,1/2)* + 

*|r 8 ,i/2) z ). 

Because of SOC, \p z Az) can be mixed with the \{p x + 
ip y ),iz) state. Similarly, the | p z ,Tz) state can be mixed 
with the | (p x —ipy)Az) state. The basis functions for Tg 
irrep, written in terms of the p orbitals, read 


1^8 5 l/^)^ d\Pz •> lA) b\ (Px T iPy)i \<z) •> 

|r s ,-l/ 2 ) z = a\p z ,i z ) +b\(p x - ipy)Az ), (B2) 


where a 2 + 2 b 2 = 1. Because the main component of 
the p orbitals is the p z orbital, we regard the parameter 
b as a perturbation. We have a b. Here we neglect 
SOC between the eg and t2g orbitals because of the large 
energy splitting between them. 

To obtain the Hamiltonian in Eq. (Bl), we transform 
the basis functions above to a set of new basis functions, 
where the spin quantization direction is along the y di¬ 
rection (the normal direction of the mirror plane). For 
each basis function in the mirror A = +i subspace, we 
have 


^(s)+i= 

1 ^ f (^0+* = ^1 dz 2 I'fy): 

V(p)+i= -F(|f 8> l/2),+i|f 8> -l/2),) = a\p z ,^ y )+ib\p x ,-\ y )+ b\p y ,iy), (B3) 


where we have used | fy) = ^(| t*) + A 4-*))- 

Thus, P and Q in Eq. (Bl) can be derived in terms 
of a, 6, and the matrix elements from k • p theory.^ The 
k • p Hamiltonian reads H k. p = — k • p = — k • (—iftV) = 

k • V, where we replace p with —ih\7 to distinguish 
it from the notation for the p orbitals. For s-p inver- 
sion, P = -i^{'H(s) +i \d z \'5’(p) +i ) = ~^a(s\d z \p z ), and 

iQ = -i£(*(s) +i \d x \*(p) +i ) = -ib£(s\d x \p x ). Recall 
that the s orbital is from Na atoms and the p orbitals are 


from Bi atoms. Similarly, for d-p inversion of NaBaBi, 
P = -i^{^(d z 2) +i \d z \^(p) +i ) = —^a(d z 2\d z \p z ), and 

iQ = - i ^('^( d z^)+i\dx\^(p)+i) = ~i^b(d z 2 \d x \p x ). 

Therefore, we have P/Q = b(sjd z !p Z ) f° r S ~P ' Ilver " 

sion, and P/Q = b()L|dx|px) ^ or i nvers i° n - This sug- 
gests that for s-p and d-p inversion, the matrix elements 

jfftfey and have °PP° site si S ns for NaBaBi. 
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